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Abstract 

We prove that the squared singular values of a fixed matrix multiplied 
with a truncation of a Haar distributed unitary matrix are distributed 
by a polynomial ensemble. This result is applied to a multiplication of 
a truncated unitary matrix with a random matrix. We show that the 
structure of polynomial ensembles and of certain Pfaffian ensembles is 
preserved. Furthermore we derive the joint singular value density of a 
product of truncated unitary matrices and its corresponding correlation 
kernel which can be written as a double contour integral. This leads to 
hard edge scaling limits that also include new finite rank perturbations 
of the Meijer G-kernels found for products of complex Ginibre random 
matrices. 


1 Introduction 

The study of random matrices benefits greatly from explicit formulas of joint 
eigenvalue densities that are known for large classes of random matrix ensembles. 
Quite a lot of these densities have the structure of a determinantal or Pfaffian 
point process. Such structures are incredibly helpful to extract the spectral 
statistics in the limit of large matrix dimension of such ensembles. Various 
techniques as free probability [54] , orthogonal polynomials [151 IHll HU ES] ; and 
supersymmetry [M] have been developed to derive these limits. 

In an important recent development it was found that explicit formulas also 
exist for the eigenvalue and singular value statistics of products of random ma¬ 
trices. This was first established for the eigenvalues [smoiiiiiHij and singular 
values 016] of products of Ginibre matrices. Shortly after this development, 
results were also derived for the eigenvalues of products of truncated unitary 
matrices [311111]. A common feature is that the joint probability densities are 
expressed in terms of Meijer G-functions which were also found in other recent 
works on random matrices mill na HIES]. 

Given the results of [311111] on eigenvalues, it seems to be natural to expect 
that also the singular values of products of truncated unitary matrices have 
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an explicit joint probability density. We aim at proving this statement. The 
squared singular values xi,... ,Xn of such a product have the joint probability 
density 

Xj) det , all Xj > 0, (1.1) 

” 3<k 

for certain functions Wk, see Corollary 12.61 below. 

A joint probability density function of the form (ED is called a polynomial 
ensemble in m- It is an example of a biorthogonal ensemble [15] whose corre¬ 
lation kernel is built out of polynomials and certain dual functions. It reduces 
to an orthogonal polynomial ensemble m in the case Wk{x) = x^ ^w{x) for 
a certain weight function w. The results of 0 E] for the singular values of 
products of complex Ginibre matrices were interpreted in |47) in the sense of a 
transformation of polynomial ensembles. Recall that a complex Ginibre matrix 
is a random matrix whose entries are independent standard complex Gaussians. 
The main result of m is the following. 

Theorem 1.1 (Theorem 2.1 in |47| I. Let n, I, V be non-negative integers with 
1 < n < 1. Let G be an {n v) x I complex Ginibre matrix, and let X be a 
random matrix of size I x n, independent of G, such that the squared singular 
values xi,... ,Xn are a polynomial ensemble 

oc - Xj) det , all Xj > 0, (1.2) 

j<k 

for certain functions defined on [0,oo). Then the squared singular 

values yi,...,yn ofY = GX are a polynomial ensemble 



j<k 

all Xj > 0, 

(1.3) 

where 



y>0. 

(1.4) 


Note that gk is the Mellin convolution of fk with the “Gamma density” x''e~^. 

We aim at an analogue of Theorem 11.11 for the product of X with a trun¬ 
cated unitary matrix and find that the structure of a polynomial ensemble is 
preserved. Instead of a Mellin convolution with a “Gamma density” we find 
a Mellin convolution with a “Beta density” a;'^(l — a;)^ with certain variables 
ly, y G No and defined on the interval [0,1], see Gorollarv 12.41 This result is an 
immediate consequence of a theorem on the transformation of squared singular 
values of a fixed matrix X when multiplied by a truncated unitary matrix that 
we present as our main theorem, see Theorem 12.II 

The spectral statistics of a generic truncation of a fixed matrix X is an old 
question. Especially in random matrix theory such truncations quite often occur 
because of its natural relation to the Jacobi ensemble m and a modification of 
the supersymmetry method m to calculate the average of an arbitrary product 
of characteristic polynomials where generic projections to lower dimensional ma¬ 
trices are needed. Applications of truncated unitary matrices in physics can be 
found in quantum mechanical evolution |33j , chaotic scattering [30j , mesoscopic 
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physics |12j and quantum information [58) . We claim that even in telecom¬ 
munications one can certainly apply truncations of unitary matrices. Usually 
Ginibre matrices model the information channel between the transmitter and 
the receiver [Si EZl EZl E]- However if the number of channels to the environ¬ 
ment is on the scale of the number of transmitting channels or even smaller 
then deviations to the Gaussian assumption should be measurable. We guess 
that those deviations result from the fact that no signal is lost. It can only be 
absorbed by the environment. Therefore the evolution in all channels should be 
unitary and, thus, in the channels between the receiver and the transmitter a 
truncated unitary matrix. 

Very recently it was shown by Ipsen and one of the authors [41] that trun¬ 
cated unitary matrices are also generally involved in the spectral statistics of 
products of rectangular random matrices, including products of matrices which 
are not Ginibre matrices or truncated unitary matrices. Moreover the trunca¬ 
tion of matrices are also important in representation theory. For example, let 
H he an m X m Hermitian matrix and T an n x m complex matrix (m > n) 
which is a truncation of a unitary matrix U G lA{m). Then the question for 
the generic eigenvalues of a truncation THT* (T* is the Hermitian conjugate 
of T) is deeply related to the fact which representations of U{n) are contained 
in a certain representation of lA{m). In particular group integrals and coset 
integrals like integrals over Stiefel manifolds are deeply related to representa¬ 
tion theory, see [33 [31 UHl I3H1 US] • Though our results are more general they 
can be partially interpreted in this framework. The relation of truncated uni¬ 
tary matrices and representation theory goes back to a parametrization of the 
unitary group by Gelfand and Zeitli rill m- They found that the eigenval¬ 
ues of THT* with n = m — 1 satisfies the so-called interlacing property, i.e. 
xi < x'l < X 2 < x '2 < ■ ■ ■ < Xm-i < 2^m-i ^ Xm where xi < ■ ■ ■ < Xm are 
the ordered eigenvalues of H and x'^ < ■ ■ ■ < x'^_-^ are the ordered eigenvalues 
of THT*. This interlacing is reflected in the “quantum numbers” labelling the 
representations of U{n) and lA{m). 

As a consequence of our generalization of Theorem 1 1.1 1 to truncated unitary 
matrices we derive the joint probability density function for the squared singular 
values of an arbitrary product Y = Tr ■ ■ ■ Ti, where each Tj, j = 1,... ,r is a 
truncation of a unitary matrix. We find a polynomial ensemble with Meijer G- 
functions which is similar to the case of products of complex Ginibre matrices [SJ 
SB]. The polynomial ensemble is a determinantal point process with a correlation 
kernel which is a double contour integral in products and ratios of Gamma 
functions and can be equivalently rewritten as a onefold integral over a product 
of two Meijer G-functions. Based on the double integral formula we are able to 
obtain hard edge scaling limits as was done in [48] for the product of complex 
Ginibre matrices. In addition to the Meijer G-kernels that are already in [IB] 
we also find certain finite rank perturbations. 

All results are summarized in section [2] The proofs of these results are 
contained in sections HHEl and ini In section |7| we briefly discuss open ques¬ 
tions to this topic. 

^We employ the transcription of Zeitlin used in m- 
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2 Statement of results 


Let us start with some preliminaries. A kxl truncation T of a matrix U G U (m) 
(m > max(fc, 1)) is a submatrix of U. Note that in the case that k + I = m the 
matrix T is equivalent with an element [T] in the coset U{m)l[U(k) x U (m — k)] 
due to the embedding 


T^[T] = \U 


- TT* 


T 


y'lm-k - T*T 


€ U{m) U €U{k) X U{m — k) 

( 2 . 1 ) 

where Ik is the k x k identity matrix, li k + l ^ m such an interpretation is not 
possible such that our studies are more general. 

A natural measure for the truncated unitary matrix T is the induced Haar 
measure of the unitary matrix U. Let dU be the normalized Haar measure. 
Then the induced measure is 


P k I 

/ T\T\S^^HTab-Uab)dU 


( 2 . 2 ) 


with + iy) = S{x)5{y) the Dirac delta function in the complex plane. 

We are interested in the singular value statistics of the matrix Y = TX where 
X G Recall that the squared singular values j/i, of the matrix Y are 

the eigenvalues of Y*Y and analogously for the matrix X. The joint probability 
density of these squared singular values involves the Vandermonde determinant 
which is given by the product A{y) = nj<fe(2/fc “ %)■ 

We present in the following subsections five results. These results comprise 
the joint probability distribution of the squared singular values of a product 
TA of a hxed matrix X (subsection 12.11) , of a random matrix X whose squared 
singular values are taken from a polynomial ensemble (subsection 12.dl) or from 
a certain Pfaffian point process (subsection 12.4L and of a product of trun¬ 
cated unitary matrices (subsection 12.Sp . Moreover we present a remarkable 
group integral involved in one of our proofs which is the analogue to the Har- 
ish Chandra/Itzykson-Zuber (HCIZ) integral [351112] for the Ginibre ensembles, 
see subsection o 


2.1 Main result 


All corollaries we present in our work are based on the following theorem. 

Theorem 2.1. Let n,m,l,v he non-negative integers with 1 < n < I < m and 
m > n-\- v Let T he an {n -\- u) x I truncation of a Haar distributed unitary 
matrix U of size mx m. Let X be a non-random matrix of size I x n with non¬ 
zero squared singular values xi,... ^Xn, and we assume them to be all pairwise 
distinct. Then the squared singular values ,..., ofY = TX have a joint 
probability density function on [0, oo)" 


oc 




det [{xk - yj)+ 


Mv) 

A(x)’ 


(2.3) 


where (x — y)+ = max(0,a; — y). The missing overall constant only depends on 
n,m and v, but is independent of X. 
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We emphasize that this theorem can be easily generalized to a matrix X 
where one or more of its squared singular values Xj coincide. Then the result 
remains valid if we replace (12.31) by an appropriate limit using L’Hopital’s rule. 

Remark 2.2. We can readily set I = n without loss of generality. The reason 
for this is exactly the same one already discussed in [5T]. We can perform a 
singular value decomposition oi X = UlX'Ub. {Ul G 1^(1), Ub. S K{n), and X' 
a rectangular matrix with only one of its main-diagonals non-zero) and absorb 
the unitary matrix Ub in T. Since I > n the matrix X' has I — n rows equal to 
zero. This structure projects the matrix T to an even smaller matrix T' of size 
(n -|- X n. Let X be the matrix X'Ub without these zero rows, in particular 
X \s n X n dimensional. Then we can consider the product Y = TX = T'X 
and apply the Theorem 12.11 for the fixed matrix X and the truncated unitary 
matrix T'. Note that X and X have the same singular values. 


2.2 An integral over the unitary group 

We give two proofs of Theorem 12.11 the first one in section H] and the second 
one in section [SJ The first proof only works in the case where m> 2n + v. In 
this case the truncation T does not have generic squared singular values which 
are equal to 1. Then T is sufficiently small compared to the dimension of the 
underlying larger unitary matrix U. In this case the set of all matrices T which 
have a squared singular value equal to 1 is a set of measure zero. 

The second proof works in all cases and is based on a calculation with test 
functions. Nonetheless, we decide to keep the first proof, too, since it is based 
on a remarkable integral over the unitary group that is of interest in its own 
right. It replaces the HCIZ integral that was used in the proof of Theorem ll.il 
The HCIZ integral formula [35l [42] is the following well-known integral over the 
unitary group U(n ), 



eyrp[tYrAUBU*]dU 



det [exp(toj6fc)]”,,^^ 
t(’^''-")/2A(a)A(6) ’ 


(2.4) 


where A and B are Hermitian matrices with pairwise distinct eigenvalues oi,..., a„, 
and bi,... ,bn, respectively, and f G C \ {0}. 

The new integral over the unitary group involves the Heaviside step function 
of a matrix argument, defined on Hermitian matrices X as 


0{X) = 



if X is positive definite, 
otherwise. 


(2.5) 


Then the generalization of (1^ for our purposes is the following theorem which 
is proven in section |31 


Theorem 2.3. Let A and B be nxn Hermitian matrices with respective eigen¬ 
values oi,..., a„ and &i,..., &„ and pairwise distinct. Let dU be the normalized 
Haar measure on the unitary group U{n). Then for every p > 0, 



det {A-UBU*Y e{A 


det 


UBU*)dU = cn,p 


(oj - bk) 


p+n—1 
+ 


J j,k=l 


A(a)A(6) 


( 2 . 6 ) 
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The constant Cn,p in (USD is 


j=0 V / 

and {uj — &fe)_|_ is {aj — bk) if the difference is positive and otherwise vanishes. 
If some of the aj and/or bj coincide we have to take the formula (12.61) in the 
limiting sense using THopital’s rule. 

Theorem l2.3l is known to hold when minofe > maxbj, see [321 formula (3.21)] 
and also [36] . Then the Heaviside step function drops out. However if this con¬ 
dition is not met it is even more surprising that the result still looks that simple, 
especially that the result can be expressed in exactly the same determinantal 
form as if the condition has been met. We emphasise that in the general case 
we usually do not effectively integrate over the whole group U (n) but a smaller 
subset. The contribution of the subset of U G U{n) for which A — UBU* is not 
positive definite vanishes due to the Heaviside step function. 

What does really happen in the integral (1^ ? To understand this we may 
choose A and B diagonal since their diagonalizing matrices can be absorbed in 
the unitary group integral. Recall that Harish Chandra [3S] traced the inte¬ 
gral (12.4p back to a sum over the Weyl group acting on B. The Weyl group 
of lA{n) is the permutation group of n elements, C U{n). This sum results 
into the determinant in (EH)- Exactly such a thing also happens here. The 
difference is that the action of w G C U{n) to the matrix M = A — ujBuj* 
sometimes yields no contribution to the integral because of the Heaviside step 
function. Because M is also diagonal the Heaviside step function of M factor¬ 
izes into Heaviside step functions of — b^(^k) for k = 1,... ,n telling us that 
Ofe — bi^(k) has to be positive definite. Despite the fact that some of the terms 
in the sum over the Weyl group vanish we can extend the sum over the whole 
group because they are zero without changing the result. This is the reason why 
inside the determinant of (12.6p we have (o^ — ^ and not (oj — bkY'^^ ^. 

Hence, one can indeed understand Theorem 12.31 bv this intuition. 

2.3 Transformation of polynomial ensemble 

Our main application of Theorem 12.II is to the situation where X is random and 
statistically independent of T, in such a way that its squared singular values are 
a polynomial ensemble on [0,oo). The proof relies on the well-known Andreief 
formula [511^. 



det [^k{.Xj)Tj^k=i 


det dp,{xi) ■ ■ ■ dp,{xn) 


= n! det 


(fk {x)/j{x)dp,{x) 


lx 


j,k=l 


( 2 . 8 ) 


that is valid for arbitrary functions ipj and Tpk on a measure space {X, p) such 
that the integrals exist. The integral (|2.8I) is used several times in our proofs. 

With the help of Andreief’s formula one can readily deduce from Theorem l2.ll 
the following Corollary. 
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Corollary 2.4. Let n,m,l,v and T be as in Theorem \2. 1\ Let X be a random 
matrix of size I x n, independent of T, such that the squared singular values 
xi,... ,Xn of X have the joint probability density function 


oc A(a;) det , all Xj > 0, (2.9) 

for certain distributions such that the moments f{x)x°‘dx with 

a = 0,... ,n — 1 exist. Then the squared singular values yi, ■ ■ ■ ,yn ofY = TX 
have the joint probability density 

oc A(y) det [gk{yj)Tj^f,^i , all y, > 0, (2.10) 

where ^ 

g,{y)= f x’'{l-xr---'-^fk(^)—, y>0, (2.11) 

Jq \X/ X 

is the Mellin convolution of fk with the “Beta distribution^^ x^{l — 
on [0,1]. 

We underline again that one can set I = n without loss of generality. 

Proof. We average (12.d|) over xi,... ,Xn with the joint probability density func¬ 
tion (El]). The n-fold integral is evaluated using (12.81) with dfj,{x) = x 
on A = [0,oo) and identifying the functions ipk{x) = fk(x) and 'f’j{x) = 
{x — ■ Then we find that the squared singular values of Y have 

a joint probability density 


oc n\ 


/ -A- \ 

■ f'OO 

n det 

/ x-^+-fk{x){x-y,)'^-^-'^-^dx 

V^=i / 

Jq 


, ( 2 . 12 ) 


J 3,k=l 


which is of the form ()2.10p with functions 

/ CX3 

x-^+^fk{x){x-yr-^-''-^dx 

Here, we applied the change of variables x i—>■ y/x. This is easily seen to reduce 
to (ifni) . □ 


Corollary 12.41 was obtained in the recent preprint [46] in a different way. 

Corollary 12.41 is the analogue of Theorem 11.11 for the case of a multiplication 
with truncated unitary matrix. It is interesting to note that Theorem 11.11 is 
obtained from Corollary 12.41 in the limit when m —^ cx) while keeping n, I, and 
r fixed, since ^JmT —>■ G where G is a complex Ginibre matrix. Recall that m 
is the dimension of the unitary matrix that T is a truncation of. 
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2.4 Transformation of PfafRan ensembles 


Theorem o can also be applied in a Pfaffian context. Instead of (12.81) we now 
use the de Bruijn formula m , see Ea Proposition 6.3.5] 



det Pf [f{xj,Xk)]2k=i ■ ■ ■ dn{xn) 


= n\ Pf 


/ / Vjix)^kiy)fix,y)d^i{x)d^i{y) 

Ux Jx 




(2.14) 


which is valid for n even and / an anti-symmetric function on X x X, i.e., 
f{y,x) = —f{x,y) for all x,y G X, such that all integrals exist. Then the 
following result is a consequence of Theorem 12.11 which can be deduced in a 
similar way as we obtained Corollary 12.41 Therefore we skip the proof and only 
state the Corollary. 

Corollary 2.5. Let n,m,l,v and T be as in Theorem \2.1\ with n even. Let X 
he independent of T such that the squared singular values of X have the joint 
probability density function 


oe A{x)Fi[f{xj,Xk)]2k=i^ all xj > 0, (2-15) 

for a certain anti-symmetric distribution f on [0, oo) x [0, oo) such that the mixed 
moments x{x 2 f{xi,x 2 )dxidx 2 exist for all j, k = 0,... ,n — 1. Then the 

squared singular values yi,. ■. ,yn ofY = TX have joint probability density 

oc A(y)Pf , all yj > 0, (2.16) 


where 


giyiiVi) 



Xl) 


m—n—v—l^v 


X2{1 - X 2 ) 


m—n — y—1 . 


\Xl X 2 J 


dxi dx2 


Xl X2 

(2.17) 


This result can be combined with Corollary 12.41 for ensembles where we 
have a mixture of orthogonal and skew-orthogonal polynomials, i.e., the joint 
probability density of the squared singular values of X is given by 


oc A(x) Pf 


f{Xj,Xk) 

-flixj) 


flixj) 

0 


1 / 

2=1,...,n 


all Xj > 0, 


(2.18) 


where n -|- n' is even and /( is an additional set of distributions. Also this 
structure carries over to the product TX as can be easily checked. 

The structure (I2.18|) is not only academically. It appears if X is real and n 
is odd (in this case we have n' = 1), see for example the real Laguerre ensemble 
or the real Jacobi ensemble E^. Also the case n' > 1 appears naturally in 
random matrix theory. For applications in QCD a random matrix model was 
proposed which breaks the Gaussian unitary ensemble by the chiral Gaussian 
unitary ensemble |20) . The joint probability density of this random matrix 
ensemble has the form (12.1811 , see El SI]■ 










2.5 Products of truncated unitary matrices 

Let us come to the product of a finite number of truncated unitary matrices. 
The joint probability density of the squared singular values readily follows from 
Corollary 12.41 As in the case of multiplying Ginibre matrices the Meijer G- 
functions play a crucial role in the spectral statistics. Meijer G-functions are 
defined via a contour integral 


G” 


Ol, 




. . . , Ltp 


\ • 1 ^m+1? • ■ • 5 


27ri 


'C 


n + s) n r(l - aj - s) 

_ 

n r(i - bj - s) n r(aj + s) 


-z-^ds, (2.19) 


where the contour G separates the poles of + s) from the poles of 

n;=ir(i — Qj — s), see musD] for more details. 

We apply Corollary [23] to the product of r truncated unitary matrices. For 
j = 1,... ,r, let Tj be a matrix of size (n + Vj) x (n + Vj-i) where t'o = 0 
and vi,... ,1'r are non-negative integers. Suppose Tj is the truncation of a Haar 
distributed unitary matrix of size rrij xrrij where mi > 2n+vi and mj > n+Vj+l 
for j = 2,..., r. The squared singular values of Ti have the joint probability 
density (|431) with parameters ni, mi and ui . This is a polynomial ensemble |25] 


1 r ci'i 1 " 

—A(a;)det wl\xY 

Zjn L 

=1 

(2.20) 

<"> = 0, 

if 0 < a: < 1, 

otherwise. 

(2.21) 


The constant normalizes the joint probability density (12.201) . Then we find 
the following Corollary. 


Corollary 2.6. Let Y = ■ ■ ■ Ti with truncated unitary matrices Tj as de¬ 

scribed above. Then the squared singular values of Y have the joint probability 
density 


UiVk-yj) det 

j<k 




1 


( 2 . 22 ) 


where w 


(r) 

k 


is given by ()2.21l) in the case r 


1 and by 



xY 




/y\ dx 
\xJ X 


(2.23) 


for j = 2,..., r when r >2. 


The weight functions are obtained as an (r — l)-fold Mellin convolution 
of the “Beta distribution”. The function can be written as a Meijer G- 
function, 


Wk\x) 


ciG 


1,0 

1,1 


( mi — 2n -|- fc 
z^i -I- fc — 1 


(2.24) 
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with Cl = r(TOi — 2n — + 1). Since the class of Meijer G-functions is closed 

under the Mellin convolution, see m formula (5)], we find from (12.231) and 


Wk\y) = 


rrir — n,..., m 2 — n,..., toi — 2n + k 

Vr,Vr-l, - ■ ■ +fc - 1 

Cr f r(^i + A: - 1 + s) nj=2 
27rz Jq r(TOi — 2n + k + s) Ylj =2 — n + s) 


y ''ds, 0 <y <1, 

(2.25) 


Here the contour G is a positively oriented curve in the complex s-plane that 
starts and ends at —oo and encircles the negative real axis. The constant Cr in 
(1^:^ is 

r 

Cr = r(mi — 2n — i^i + 1) r(TOj — n — Vj). (2.26) 

i=2 

It can be checked from ()2.25l) that the linear span of ,..., Wn'^ consists 
of all functions of the form 


1 f g(5)nj^ir(s + i^j) _ 
^ni Jc + mj - n)^ 


0 < y < 1, 


(2.27) 


where q{s) is a polynomial of degree smaller than n. Remarkably enough, this 
space does not depend on the ordering of the parameters mi,, mr and neither 
on the ordering of the parameters vi,... ,Vr. 

The fc-point correlation function of the joint probability density (12.221) satis¬ 
fies a determinantal point process on [0,1]. In the center of such a determinantal 
point process stands a correlation kernel that can always be written as 


n—1 

Kn{x,y) = ^ Pk{x)Qk{y). 

k=0 


(2.28) 


For a polynomial ensemble, the function Pk is a polynomial of degree k and Qk 

(r) (r) 

is in the linear span oi w\ ,..., Wn satisfying the biorthogonality 


/ Pj{x)Qkix)dx = Sj^k, 

Jo 


(2.29) 


see e.g. m- As in [IS] we find integral representations for the biorthogonal 
functions Pk and Qk and a double integral formula for Kn. In what follows we 
use the Pochhammer symbol 


{a)k = a{a -I- 1) • • • (a -I- fc — I) = 


r(a -I- k) 


r(a) ■ 

Proposition 2.7. For every k = 0,... ,n — 1, we have 


(2.30) 




1 


r(i + 1 + TTlj — Tl) 


_ ^0,r+l 

~ ^r+l,r+l 


{t-k)k+i^^^ T{t + l + Vj) 

k + l,n — mi,... ,n — mr 

0,-vi,... ,-iyr 


* dt 


(2.31) 
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where E/c is a closed contour encircling the interval [0,fc] once in positive di¬ 
rection and is not encircling any pole of the integrand in (12.311) that is outside 
[0,fc]. Moreover we have 


Qk{y) — 



r(s + nij — n) 


y ‘'ds, 


^r+1,0 

^r+l,r+l 


—/c, mi — n,..., TUr — n 


0 < 2 / < 1 , 

elsewhere 


(2.32) 


with the same contour C as in (12.251) . 
The kernel Kn from (12.281) is 


Kn{x,y) 


1 

(27ri)2 


Jc r(t + 1 + u,)r{s + l + m,-n) s-t 


— I ^0,r+l 

/ '-^r+l,r+l 


n — mo, . . . , 71 — TUr 


UX 


^ '-^r+l,r+l 


mo — n,..., rur — n 


Z/Q,.- 




uy I du 


(2.33) 


which is valid if E^, and C do not intersect. In (12.331) it is understood that 
mo = Vo = 0 . 


The proposition is proved in section [51 


2.6 Hard edge scaling limits 

Based on the double integral representation (I2.33P we analyze some scaling limits 
of the correlation kernel as n —>■ oo. In a forthcoming publication we show that 
the usual sine kernel limit can be found in the bulk and the Airy kernel limit at 
the soft edges. These results will be reported elsewhere, see also [49] . 

Here we want to look at a more straightforward scaling limit, namely the 
hard edge scaling at the origin in the following situation. Taking n —>■ oo, we 
simultaneously have to let mj —>■ oo for j = 1,..., r, since mi > 2n vi and 
mj > n Vj 1 ioT j > 2. We keep Vj fixed for every j = 1,... ,r. We choose 
a subset J of indices 

J = C {2,...,r}, with 0 < 9 = |J| < r (2.34) 

and integers pLi,..., y,q with > Vj,, + 1, and we assume 

mj — n —>■ oo for j G {1,..., r} \ J, 

mjj^ — n = pik \s constant for j = jk € J. 

This leads to our final result, which we also prove in section |5| 

Theorem 2.8. In the above setting, we put 

Cn = nY^imj — n). (2.35) 

iiJ 
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Then the kernels Kn from (I2.33|) have the following hard edge scaling limit, 


hm l/fji JL 

n->oo c„ \ C„ C„ 

_ 1 ^ 2+®°° I" ___ + 1 + Vj) sinTTs -pr r(i + 1 + /ife) x*y ^ 

ifl-nif ^ is r(i + 1 + sinTTt T{s + l + fj,k) s-t 



^g,r+l 


/il, . . . , /iqr 





Ml; ■ • ■ 5 Mg 


(2.36) 


where H is a contour around the positive real axis in the half-plane Re i > — ^. 

The kernel (12.361) reduces to the Meijer G-kernel described in [?5] in case 
q = 0. These kernels appeared before for limits of products of Ginibre matrices 
[48) , products with inverses of Ginibre matrices [26] , for biorthogonal ensembles 
m, and also as limits for Cauchy two matrix models [HillllH]- According 
to Theorem 12.81 we obtain the same limits for products of truncated unitary 
matrices provided that the dimensions ruj of the underlying unitary matrices 
become large compared to n, in the sense that mj — n —>■ +oo for every j. 

The kernels (12.3611 are new for q > 1, and these are finite rank perturbation 
of the Meijer G-kernels from [45]. To see this, we recall that f^k > + 1. Let 

us assume, for notational simplicity, that jk = r — q -\- k in (12.341) . Then 


Rit) n 

/c=l 


r(t +1 + /ifc) 

r(t + 1 + Ur-q+k) 


(2.37) 


is a polynomial of degree deg R = — Vr-q-i-k), and (12.3611 can be written 

as 


2 r 2 +'^^^ r ^{s1Vj) siniTs Rft) x*y ® ^ 
(27ri)2 Mi-ioo r(t+1 + ^'j) sin TTt i?(s) s-t 


(2.38) 


This is indeed a finite rank perturbation of the Meijer G-kernel with parameters 
vi,, Vr-q, since i? is a polynomial. In particular for g = r — 1, it is a finite 
rank modification of the hard edge Bessel kernel. Such finite rank modifications 
were also obtained in [23] in a somewhat different context. 

In the ensuing sections we prove our statements. We start in section jSj with 
the proof of Theorem 12.31 since it is used in the first proof of Theorem 12.11 that 
we present in section |4| The second proof is in section |5| This proof is a rather 
lenghty sequence of matrix integral evaluations and we have broken it up into 
six steps. The proofs of Proposition l2.7l and Theorem l2.8l are shown in sectionlHl 


3 Proof of Theorem 12.3 

For the proof of Theorem 12.31 we need the Ingham-Siegel formula [351 US] 
I expliTiHX]det{H — zIn)~^~^dH = cexpjfzTrX] detX^6((A) 


«(») 


(3.1) 
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with a normalization constant c depending on p and n, only, which can be fixed 
by the choice X = In and z = i. This integral is valid for Hermitian matrices 
X € 'H(n) and Imz > 0. The integral is over the space 'H(n) = g[(n)/u(n) of 
n X n Hermitian matrices H with the flat Lebesgue measure, 

n 

dH = n dHjj n dUeHjkdlmHjk- (3.2) 

j=i j<fc 


Here 0 ((n) and u(n) are the Lie algebras of the general linear and the unitary 
group, respectively. If p is not an integer then we define via the spectral repre¬ 
sentation, 

n 

det{H - zIn)-^-P = l[{h, - z)-"-P, (3.3) 

i=i 

where hi,... ,hn are the real eigenvalues of the Hermitian matrix H and {h — 
is defined in the complex h-plane with a branch cut along {z+iy | y > 0 }, 
and with h”+P(/i — ^ 1 as h —^ -l-oo. 

For n = 1 the Ingham-Siegel formula (ED reduces to 



^ixs 

7- zn—ds = 

(s - zy+p 


2TTie^^P/^ 


r(p- 

.0, 


1 ) 


if a; > 0 , 
if a; < 0 , 


(3.4) 


which is valid for Imz > 0 and p > 0. It can be verified by contour integration. 
Proof of Theorem 12.31 We consider the integral 


J{A,B)= det{A-UBU*f0{A-UBU*)dU. (3.5) 

Ju{n) 

The determinant can be rewritten via the Ingham-Siegel formula (ED identify¬ 
ing X = H - UBU*. We obtain 

J(H,H) oc / f y'^’^^^^-^^^^'>det{H-zIn)~^~PdHdU 

Ju(n) J'H{n) 

(3.6) 

Both integrals are absolutely integrable because the integral over U is over a 
compact set with a continuous integrand and the integral over H is bounded by 
|A(/i) nj=i(^i “ z)~^~P/A{a)\ for o®- Recall that oi,..., a„ are 

the real, pairwise distinct eigenvalues of A. Hence we can interchange the order 
of integration. 

The integral over the unitary group is evaluated with the HCIZ formula 
(12.41) . Then (|3.6I) is up to a constant 

detiH - (3.7) 

Juin) A{h)A{b) 

We write H = VhV* for the eigenvalue decomposition of H where V G U{n) 
and h = diag(/ii,..., /i„). Moreover we use that dH oc A(h)^ dVdhi... dhn, see 
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e.g. m- Using this in (13.71) leads to 


J{A, B) oc 


^—iz Tr{A — B) 


iTrVhV*A 




—n—p 


JR" JU(n 

X A{h) det dVdhi ■ ■ ■ dh. 

The integral over V G 14(n) is again a HCIZ integral (j2.4p , 


(3.8) 


J{A,B) oc 


„-iz Tr(A-B) 


det [( 


^ihj flfe 1 


A(a)A(6) l 

xdet 


i=i 
dh^i . 


—n—p 


(3.9) 


The factors in the product 0^=1 — z) ^ ^ can be pulled into either one of 

the two determinants. Then, Andreief’s identity (12.81) can be applied to find 


J{A,B) oc 


g —22 Tt(A — B) 

A(a)A(6) 


det 



^is(aj-bk) 

(s - z)"+P 





(3.10) 


The integral in the determinant is of the form dH which is up to a constant 
equal to — 6^)"^^”^. The exponential factors e'^Aa-j-bk) inside the 

determinant cancel with those in front of the determinant. The resulting ex¬ 
pression is 


J[A,B) oc 


det 


bk) 


n+p—1 
+ 


A(a)A(6) 


(3.11) 


which is up to a constant indeed the right hand side of (12.61) . 

Whenever A — UBU* is positive definite for all U S U{n), the formula (12.61) 
reduces to 



det 

det (A - UBU*f dU = Cn,p - 


A(a)A(6) 


(3.12) 


This is equivalent to an integral given by Gross and Richards in |321 formula 
(3.21)], namely 


det 


(1 Sjtf,) 


j,k=l 


A(s)A(t) 



517^7*)"^“+”"^^ dU 


(3.13) 


with Cn,a = nj=o^(®)i/j'- The formula (I3.13P is valid whenever S and T are 
Hermitian matrices with eigenvalues si,..., s„ and ti,... ,tn, respectively, sat¬ 
isfying \sjtk\ < 1 for all 7, fc = l,...,n. The formulas (13.121) and (13.131) are 
related if we take —a = p + n — 1, S = A~^ and T = B. The constants are re¬ 
lated by Cn,p = ilcn,a, and the formula (|2.7I) follows which completes the proof 
of Theorem l2.3l □ 
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4 First proof of Theorem 12.1 

As already said, our first proof of Theorem 12.11 only works if m > 2n + In 
that case there is an explicit formula for the distribution of a truncation T of 
size {n + v) x n, namely 

c det(/„ - - T*T) dT (4.1) 

with dT = d Re Tjk d Im Tjk the flat Lebesgue measure on the space 

of (n + zz) X n rectangular complex matrices and c is a constant. The function 
9 is the Heaviside step function of a matrix argument defined in ()2.5I) . The 
(unordered) eigenvalues ti,..., of T*T are in the interval [0,1] and have the 
joint probability density 

n 

0<ti,...,t„<l. (4.2) 

^, -1 

j<k 3 = 1 

This is an example of a Jacobi ensemble [18]. When m <2n + iy, then T*T has 
the eigenvalue 1 with multiplicity of at least 2n + z/ — m and the density dSD 
is not valid anymore. 

We follow the proof of Lemma 2.2 in m except that at a certain stage in 
the proof the HCIZ integral (|2.4I1 is replaced by the integral (12.611 . 

Proof of Theorem 12.11 in the case m>2n + v As already discussed in Re¬ 
mark [52] we may restrict to the case I = n without loss of generality. 

Consider a fixed square matrix X of size n x n which is assumed to be 
invertible. The change of variables T = TX has the Jacobian, see e.g. m 
Theorem 3.2], 

n 

det(A*A)-("+") = Y[ (4.3) 

k^l 

The distribution (SU) on T (where T has size (rz -|- zz) x n) then transforms into 
the distribution 

n 

oc Y[ det(/„ - {X-^)*Y*YX-^)^-^^-''e{Ir^ - {X-^)*Y*YX-^) dY 

k^l 

71 

= det(X*A - Y*Y)^-'^^-''e{X*X - Y*Y) dY. (4.4) 

k=l 

In the next step we perform a singular value decomposition Y = VYU with 
Jacobian [24] 

n 2^1 j MyfdUdVdyi ...dy^ (4.5) 

written in terms of the squared singular values yi,... ,yn of T. The measure 
dU is the Haar measure on lA (n) and dV is the invariant measure onU{n + i')/ 
[W"(l) X We use this fact in (14.411 to perform the integration of V, 


dY ^ 
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which only contributes to the constant. This yields a probability measure on 
U{n) X [OjOo)" proportional to 


oc 




det(X*X - 

X e{X*X - U*Y?U)X{yf dUdyi ■ ■ ■ dy^. 


(4.6) 


The integral over U in (ITHI) can be done with the help of the integral (ESI) 
with A = X*X, B — Y?, and p = m — 2n — v. This leads to the density for the 
squared singular values yi,... ,yn of 1", given those of X, which is proportional 
to (12.3|) and Theorem ED follows for m>2n + v. □ 


5 Second proof of Theorem 12.1 

We underline that the second approach to prove Theorem 12.11 does not rely 
on the restriction m > 2n + v. As before we denote the set of n x n unitary 
matrices and of n x n Hermitian matrices by U(n) and 'H{n), respectively. We 
also use M{m,n) for the space of m x n complex matrices and abbreviate 
M{m) = M{m,m). 

Also in the second approach we assume that I = n because it does not 
restrict generality, see Remark 12.21 We assume A to be a fixed n x n matrix 
with non-zero squared singular values. 

5.1 Preliminaries 

Let / be a symmetric function in n variables. We extend / to Hermitian matrices 
A by defining f{A) = /(oi,..., a„) if oi,..., a„ are the eigenvalues of A. With 
a slight abuse of notation we also define f{B) for (n + x {n + v) matrices 
B having v eigenvalues equal to 0, by putting f{B) = /(6i,..., bn) if 6i,..., 
are the non-zero eigenvalues of B. 

Then our aim is to prove that for all continuous symmetric functions / on 
[0, oo)", we have 

E[f{Y*Y)]= f fiyi,...,yn)px,Yix,y)dyi...dyn (5.1) 

J [0,oo)^ 

where, for a given X, y ^ px,Y{x,y) denotes the density from (12.31) . It will be 
enough to prove ED for symmetric polynomial functions /, since the density 
Px,Y has a compact support and the symmetric polynomials are then obviously 
uniformly dense in the set of all continuous symmetric functions. 

Note that, by our definition of / on matrices, we have f(Y*Y) = f(YY*) = 
f {TXX*T*). Since T is the truncation of a Haar distributed unitary matrix U 
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of size m X we have 


IU{m) 


E[f{YY*)] = / (/„+, 0„+,, UXX*U* 


^n-\-u 


Ora—n — u,n-\-h> 


du 


[ f( (/„+, On+.,m-u-.) MXX*M* 

J Ai{m,n) \ 




X Yl d<^^\{M*M}jk)Y[6{{M*M},,-l)dM, (5.2) 

l<j<k<n j—1 


where dU is the normalized Haar measure on the unitary group lA(rn). The 
matrix Op^q is the zero matrix of size p x q. The complex m x n matrix 


U = U 



(5.3) 


is an element in the Stiefel manifold lA{m)/lA{m — n). The orthonormality of 
the columns of U can be enforced by Dirac delta functions (recall that is 
the one for complex numbers). In this way we integrate in (ESI) over the larger 
space M{m,n). See also the discussion in [TO] . 

The complex matrix M can be partitioned into two blocks 

with Ml an (n + v) X n complex matrix and M 2 an (m — n — v) x n complex 
matrix. Then we have to calculate 


E[f{YY*)] 


/ A4{n-\-iy,m) 


f {MiXX*M*) 


n+i>' 


X n S^^K{M*M}qk)l[S{{M*M}, 
l<j<k<n-\-i' j—1 


l)dM. 

(5.5) 


5.2 Proof of Theorem 12.11 

To establish (ED we proceed in six steps. 

Step 1: Matrix delta function In the first step we rewrite the Dirac delta 
functions in (ED as Fourier-Laplace transforms m where 

n 

n ({ M * M }, fc ) n k{m*m},, - 1 ) 

l<j<k<n j—1 

= lim [ exp[Tr(/„ - zif)(/„ - M*M) -tTvK^jdK (5.6) 

t^o2"7r"- 

with an integration over the space 'H{n) of Hermitian n x n matrices K. For 
an integration over the whole group U{m), i.e. m = n, this integration is equal 
to the one in EZl formula (13)]. The shift of the matrix K by ensures the 
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absolute integrability of the integral over M and the Gaussian incorporating the 
limit in the auxiliary variable t guarantees the absolute integrability in K. Note 
that the limit has to be understood as a limit in the weak topology meaning 
that we have to integrate over M, first, and, then, take the limit t ^ 0. Hence 
the integral dES]) reads 

E[f{YY*)](x\im [ f{MiXX*Ml) 

JM{m,n) 

X / exp[Tr(/„ - iK){In - M^Mi - M 2 *M 2 ) - t Tr K^]dKdM (5.7) 

up to a constant only depending on m, n and ly. 

Both integrals are absolutely integrable. Therefore we can interchange the 
integrals. The integral over the matrix M 2 is a Gaussian integral yielding 


E [f{YY*)] oclim / [ f (MiXX*M*) det(/„ - iK)-^+'^+'' 

J-Hin) JM{n+v,n) 

X exp[Tr(/„ - iK){In - M*Mi) - t Tr K^]dMidK. (5.8) 

Finally we take t —>■ 0. This can be done because / is polynomial. In the case 
/ = 1 the integral over Mi yields an additional factor det{In+u — iK)~'^ ensuring 
the absolute integrability also at t = 0. Since the function / is polynomial we 
know that f [MiXX*Ml) is a polynomial in the matrix entries of Mi. Thus 
the Gaussian integral over Mi yields a polynomial in {In+v — iK)~^ with the 
lowest order to be det(/„+,y — Therefore the integrand of the iF-integral 

after integrating over Mi, first, is indeed absolutely integrable also at t = 0 such 
that 

E[f{YY*)](x [ [ f{MiXX*Ml)det{In-iK)-”^+'^+'' 

X exp[Tr(/„ - iK){In - MlMi)\dMidK. (5.9) 

We underline that now the order of the integrals is crucial and cannot be inter¬ 
changed. 

Step 2: Changes of variable The change of variables Mi 1 —>■ S 
the Jacobian 

n 

det{X*X)-^-'' = 

and the change of variables K 1 —>■ K = X~"^K{X~^)* on the space of Hermitian 
matrices yields 

71 

det{X*X)^ = Y[x]. (5.11) 


= MiX has 
(5.10) 
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Applying this to (15.91) (and dropping the tilde from K), we obtain 

n « 

E [f{YY*)] (xYlx-'' dKdet{In - 

Z-I Jnin) 


f=l 

X 




dSf{SS*)e' 


n « 

TT / dKdet{{X*X)-'^ - 

Jwn) 


f=i 

X 


/ Ai{n+i',n) 


lH(n) 


j=l JH{n 

where dt is defined by 
d!{K) := det(/„ - 


n p 

TT^-m+n/ ^^J(K + ^{X*X)-^-iIn)dK, (5.12) 

JWn) 




(5.13) 


Step 3: Shift in the matrix K In this step we prove the following Lemma. 

Lemma 5.1. For any matrix A with ImA := [A — A*)/2i > —(which means 
that In +ImA is a positive definite Hermitian matrix), we have 


I d>{K + A)dK = [ d>{K)dK. 


(5.14) 


Proof. If A is Hermitian then we can simply apply the linear change of variables 
AT + A i—7> AT to obtain (|5.14|1 . Therefore we may restrict to the case A = iB 
with B Hermitian and H + /„ positive definite. We may also restrict to the case 
where H is a diagonal matrix. To see this we write 'Lx to indicate that the 
definition (|5.13l) depends on X. Then for a unitary matrix U one has 


d>x{UKU*) = 'i>u.x{K), (5.15) 

which follows from inserting UKU* into the definition (15.131) and changing vari¬ 
ables S !->■ SU*. Recall that f{SS*) = f{USS*U*) for all U € U{n) because / 
only depends on the squared singular values of S. The invariance (15.151) implies 
by the unitary invariance of dK, see e.g. [21], that 



d!x{K + iB)dK 



d!u-x{K + iU*BU)dK. 


(5.16) 


and we may choose the unitary matrix U so that U*BU is diagonal. 

Let p = m — n — v>Q and A = iB = i diag(6i,..., bn) a diagonal matrix 
with bj > —1, for j = 1,..., n. Note that for 0 < t < 1, 

|4'(A: -f itB)\ < |det(/„ +tB- iK)\~P [ \fiSS*)\ eT'Hn+‘S)(xx*-s*s)^ 5 . 

Juin) 

<Co\det{In + tB-iK)\~^. ( 5 . 17 ) 
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The prefactor is a finite constant Co > 0 depending on B and X but is inde¬ 
pendent of t G [0, !]■ Since all hj > —1, it is clear that /„ + tB is a, positive 
definite matrix. Hence there exists a constant Ci > 0, independent of t, such 
that |det(/„ -I- tB)\ > Ci resulting in 


\det{In + tB - iK)\ > Cl |det(/„ - iMt)\ , (5.18) 

where 

Mt:=iIn+tB)-^K{I„+tB)-^ (5.19) 

is a Hermitian matrix. Let Xi{t),, Xn{t) denote the eigenvalues of Mt. Then 
it is easy to check that, since the eigenvalues Xj{t) are real, 

n n 

|det(In - iMt)\ = n = n + 

k^l k^l 

> max |Afe(t)| = max ^)l (5.20) 

fc=l,....n a;GC"\{0} {x,x) 

by the properties of Rayleigh quotients. Taking x = xj = {In +tB)i ej in (15.2011 
and noting that {MtXj,Xj) = Kjj, see (15.1911 . we obtain 

|det(/„ — jMt)| > max ^ (5.21) 

J = l,...,n {Xj,Xj} 

The vectors Xj depend on t, but their norms are uniformly bounded from below 
(xj , Xj) >1-1- t6min > C > 0 for t G [0,1] where 6niin = minj=i_..._„ hj and C = 1 
if ^min > 0 or C = 1 — 6niin if 0 > > — 1- Then, combining (I5.17L (15.1811 

and ()5.21L we obtain 


\^{K + UB)\ < 


C 2 

maxj^i,„„„ \Kjj\p' 


(5.22) 


for some constant C 2 > 0, independent of t. 

In the integral 'd>{K)dK we first do the integration over the diagonal 

elements Kjj for j = 1,... ,n. The integrand is analytic in each of the Kjj. 
Because of the estimate ()5.22|l with p > 1, the integral can be deformed from 
the real line to the horizontal line in the complex Kjj-plane with imaginary part 
equal to bj. We do this for all diagonal entries resulting in (15.1411 . □ 


Lemma o can be applied to (I5.12p because A = i{X*X) ^ — iln satisfies 
luiA = {X*X)~^ — In > —In- Hence we have 


n p 

E[f{YY*)] oc n / dK det{In - 

j=i Jnc 


f A1(n+i/,n) 


dSf{SS*)e 


*\„TT{Ir,-iK)(X‘‘X-S*S) 


(5.23) 


Step 4: Singular value decomposition of S We take the singular value 
decomposition S = V 1 EV 2 where Vi G U{n + i/)/[W”(l) x U{iy)], V 2 G U{n), 
and E is a diagonal matrix with the singular values of S. The Jacobian of this 
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transformation is proportional to A(?/)^ 0^=1 V’j^ see [23], where yi,. ■ ■ ,yn are 
the squared singular values. Recall that / is defined as a symmetric function 
on the eigenvalues. Then (15.131) reads 


'S{K) oc 


gTr(7„-iiC)(X*X) 

det(/„ - 


« n 

/ fiyi,...,yn)T\yjMyfdyi...dyn 

ifO.oo)" 

X [ (5.24) 


'lA{n) 

The integral over Vi^U (n) is a HCIZ integral (12.4|) . Let Ai,..., A„ be the eigen¬ 
values of RT, then the HCIZ integral yields a term proportional to det[e“^^“*''*-’^^'‘] 
/[A(A)A(j/)]. We end up with 


4'(R:) 


pTr(/„-i7f)(X*X) 


det(/„ - iR:)™-"-‘'A(A) 

P n 

/ f{yi,---,yn)Wy''Hy)det 

'[O.oo)- 


= -(l-*^3)2/fc 


J j,k=l 


dyi...dyn. (5.25) 


Step 5: Eigenvalue decomposition of K Next, we decompose K = VqAVq 
in a unitary matrix Vq G K{n)ll{"'{l) and its eigenvalues A = diag(Ai,..., A„). 
The Jacobian is proportional to A(A)^, see e.g. |2I], and by (15.231) and (15.2511 
we have 



The Vo integral is again a HCIZ integral (Elj), and it gives a contribution propor¬ 
tional to det[e(^“®'^J^“^'=]/[A(A)A(a;)]. Plugging this term into (15.2(111 and noting 
that we may change the order of integration at this stage, we find 


E[f{YY*)]= f{yi,---,yn)px,Y{x,y)dyi---dyn (5.27) 

J [0,oo)’^ 


where 


Px,Y{x,y) oc 




A(a:) 






n 

det 


- 

jA=i 

■ 


dAi•• • dXn- 


(5.28) 
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Step 6: Andreief formula We calculate the integral over Ai,..., A„ in (15.281) 
via the Andreief formula (12.8p . This gives the determinant 


nl det 



(1 - 



3,k=l 


(5.29) 


with integrals of the form (13.41) . The result is 



75-T5- dX 

(1 - zA)™-"-*" 


27r 

(to — n — — 1)! 


[xk -yj)+ 


m—n—i^ — 1 


(5.30) 


We use this in (15.271) and (15.281) and obtain (15.11) . This concludes concludes the 
proof of the Theorem 12.II □ 


6 Proofs of Proposition 12.71 and Theorem 12.8 

We first prove Proposition 12.71 and Theorem 12.81 afterwards. 


Proof of Proposition 12.71 Recall that Pk and Qk are defined by (12.311) and 
(I2.32|) . respectively, for k = 0, ...,n — 1. Then Qk is in the linear span of 
w^\ ..., Wn'^ since it is of the form ()2.27l) with q{s) = (s — k)k a polynomial of 
degree k < n—1. Also note that the integrand in ()2.3ip has poles at t = 0,..., fc, 
and no other poles inside E^. Thus if we evaluate (12.311) by means of the residue 
theorem and find contributions from t = 0,..., fc, only. Hence the function Pk 
is a polynomial of degree k. It remains to verify the biorthogonality (12.291) . 

Let I, k = 0,1,... ,n — 1. From (12.311) we have 


Pi{x)Qk{x)dx 


n 


r(t + 1 + rrij — n) 


Isi (t ~ Oi+i r(t + 1 + Vj) 


x*'Qk{x)dxdt. (6.1) 


The moments of Qk are given by the general identity for the moments of the 
Meijer G-function, 


.1 r 

/ x‘^~^Qk{x)dx = {s - k)k\\ 
do 


r(s + Vj) 


r(s + TO,- — n) 

This identity can be plugged into dSU which cancels a lot of T-factors, 

r[t-z] 


( 6 . 2 ) 


Pi{x)Qk{x)dx = ^ j) 


it + l-k)k^^^ 1 


(t - l)i+i 


27rz r[t + 1 — A:] 


dt. (6.3) 


For I < k, the integrand in (lOl) is a polynomial and the integral is zero by 
Cauchy’s theorem. For I > k, there are poles at t = k,... ,l which are all inside 
the contour E/. The integrand is a rational function that behaves like 
as |t| —^ oo. Therefore by simple residue calculation at infinity, the integral also 
vanishes if 1 > fc and it is equal to 1 if / = fc. Thus (12.291) is satisfied. 
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Inserting (12.311) and (12.321) into (12.281) , and noting that we can take the same 
contour for every k in (12.311) . we obtain a double integral for Kn, 


Kn{x,y) = 

1 


1 


(27rz)^ 


(27ri)2 

ds 


[ ds I dt'n r(s + ^'j-)r(f+ l + mj- -n) (s - k)k -s 

Ic + -n)T{t+l + iyj) k)k+i 

TT r(s + 1 + Uj)T{t + I + ruj - n) (s + 1 - k)k t -s-i 

J-A T/'c-Ul -L m ■ — -L 1 /A (i — 


fC 


7=1 r(s + 1 + TTij - n)r(t + 1 + Uj) ^ (i - k)k+i 


(6.4) 


where we made the change of variables s i—>■ s + 1. The summation can be 
simplified, because of the telescoping sum 

•^ {s + l — k)k _ 1 / r(s + 1) r(f + 1 — n) \ 

^ "s-tVr(t + i)r(s + i-n) 


which can be readily checked by complete induction. The contours and C 
do not intersect. Therefore the term t = s is not a pole inside the contour I]„. 
The double integral (16.4p splits into two terms due to (j6.5l) . For the second term 
the contour integral over t does not encircle any pole because of mj — n > Vj 
and rrij — n, Vj £ No- Thus this term vanishes. The remaining term is the one 
shown in (I2.33|) . 

The expression (12.331) in terms of Meijer G-functions is obtained by noticing 
that 

-= - / u^-^-^du. ( 6 . 6 ) 

s — t Jq 

Interchanging the s and t integral with the u integral and using the definition of 
the Meijer G-functions (I2.19P the identity follows, see also the proof of Theorem 
5.3 in [45]. This concludes the proof. □ 

Proof of Theorem 12.81 We employ the following asymptotic behavior of a ra¬ 
tio of Gamma functions 


r(t -I- 1 -I- rrij — n) 
r(s + 1 + mj — n) 


sinTTS 
sin TTt 

{rrij - 


(^1 + O ( 
n)‘-" (l + O 



for j = 0, 


for j G {l,...,r}\ J 
(6.7) 


as n —>■ 00. This follows as in the proof of Theorem 5.3 in [45], since toq = 0 
and rrij —n—j-ooasn—?>oo for j G {1,..., r} \ J. 

In the double integral formula in (|2.33l) we deform the contour to S, and 
obtain 


—Kn 

Cn 


Cn Cn J 


(27rz)2 



r(s + -p|- 

r(t -I-1 -I- Uj) 


r (t -I- 1 -I- rrij 
r(s + 1 -I- rrij 


n) x^y ® ^ 
n) s — t 


( 6 . 8 ) 
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Because of definition (12.351) we have 


s-t TT r(f + \ + mj — n) 

" f},r{s + l + m,-n) 


'^^-t r(t + l + TOo-n) \ 


r(s + 1 + mo - n)/ 


X 



r(s + 



Each of the factors in the product has a finite limit as n —>■ oo, cf. dlTTl) . and 
the full product tends to 



( 6 . 10 ) 


sinTrt ^^r(s + l + /rfc) 


since mj — n = /Xfc if _) = jk G J- It is then allowed to take the limit inside the 
double integral in (16.81) and deform the contour C to the vertical line Res = 
— 1/2 by analyticity and the decay of the integrand at infinity in the s-variable. 
This leads to (|2.36l) . 

The expression in terms of Meijer G-functions is obtained by using the re¬ 
flection formula of the Gamma function, sinTrz = 7r/[r(l — z)T{z)] together 
with (EH). The identity now follows along the same lines as in the proof of 
Proposition 1 2. 71 □ 

7 Conclusions and Outlook 

We analyzed the singular value statistics of a matrix X (fixed as well as ran¬ 
domly chosen) multiplied by a truncated unitary matrix T which is distributed 
by the induced Haar measure. Though we only considered a multiplication from 
the left side TX one can easily generalize the results to a multiplication from 
both sides TlXTr where Tl and Tr, are truncations of two independent unitary 
matrices. The reason for such a simple generalization is the determinantal point 
process fulfilled by the joint probability density of the singular values. We proved 
that the joint probability density of the squared singular values of TX satishes 
a polynomial ensemble if the joint probability density of the squared singular 
values of X does. In particular with the help of our results one can calculate 
the squared singular value statistics of any product Tl,i • • • TL^rL-^^R.i ‘ ’ ’ TR,rR 
with rLjTR S No and Tlj and Trj truncations of independent unitary matrices 
and X either fixed or another random matrix. Indeed one can also mix prod¬ 
ucts of truncated unitary matrices with Ginibre matrices. The combination of 
Theorem o and Corollary 12.41 yields a polynomial ensemble for the squared 
singular values of such a mixed product. We expect that also for such a mixed 
product the statistics are governed by Meijer G-functions since this particular 
kind of functions is closed under Mellin convolution as it was shown here for a 
pure product of truncated unitary matrices and studied in for a product 

of Ginibre matrices. 

Our study shows that the determinantal point process also applies to a prod¬ 
uct with truncated unitary matrices. In particular one needs a group integral 
replacing the Harish Chandra/ Itzykson-Zuber integral [35l |42] for which it is 
not immediate that the result can be written in terms of determinants, cf. The¬ 
orem o This particular result is even more astounding when noticing that 
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one does not effectively integrate over the whole unitary group but only over a 
subset. The reason for this is a Heaviside step function in the integrand. Har- 
ish Chandra made contact between group integrals and representation theory 
m- It would be interesting if also something like this exists for the integral 
considered by us and can be explained by group theory. 

Additionally we looked at the spectral statistics of a product of truncated 
unitary matrices in detail. For this product we calculated the kernel of the cor¬ 
responding polynomial ensemble of the squared singular values at finite matrix 
size and in the hard edge limit at the origin. In a forthcoming publication we 
also derive the bulk, the soft and the hard edge statistics as it was very recently 
done for the Ginibre ensemble in [49]. The latter two statistics may appear 
at the lower and the upper bound of the support because the squared singular 
values of a product of truncated unitary matrices always live on the interval 
[0,1]. If the support touches either the origin or the upper bound 1 one would 
expect hard edge statistics at these points. 

Another generalization of our results refers to the restriction that the first 
matrix in the product of truncated unitary matrices has to satisfy mi > 2n + iyi, 
see Corollary 12.41 This matrix has not generally to be the first matrix Ti. 
With the help of the discussion above it can be any matrix in the product. 
Nevertheless we have to assume that at least one truncated unitary matrix 
multiplied has to satisfy the condition to prove Corollary 12.41 in the way we 
have done. An interesting question would be: What happens if this restriction 
is not met? From the spectral statistics of one truncated unitary matrix, see e.g. 

we know that some singular values are exactly located at 1. Numerical 
simulations performed by us hint that this seems to be true also for a product 
of truncated unitary matrices. In the case that this is indeed true, the question 
arises about the algebraic structure. Does the determinantal point process carry 
over to a product ■ ■ - Ti where the restriction mj > 2nj + Vj is not met for 
all j = l,...,r? 
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